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Description 



Retrospective andAnd/Oor Wwindow-Ceontrolled Ffiltering Oof Kmages Ffor 
adopting A dapting Ss harpness Aftnd Nnoise I*n Ceomputetd Ttomography 

fOOOn Th e present application hereby claims priority under 35 U.S.C. 8119 on 
German pateni application number DE 10238322.7 filed August 21. 2002. the 
entire contents of which are hereby incorporated herein bv reference. 

BACKGROUND OF THE INVENTION 

100021 Image data of a measurement object to be examined can be obtained with 
the aid of modern medical diagnostic methods such as, for example, X-ray 
computer computed tomography (CT). As a rule, the measurement object 
examined is a patient. 

£0.003] X-ray oomputo r cojnputed tomography - noted below by CT for short - is a 
specific X-ray recording method thai differs fundamentally in image structure from 
the classical X-ray tomography method. What is obtained in the case of CT 
photographs are transverse sectional images, that is to say images of body slices 
that are aligned substantially perpendicular to the body axis. The tissue-specific 
physical quantity represented in the image is the distribution of the attenuation 
value of X-ray radiation m(x, y) in the sectional plane. The CT image is obtained 
by reconstructing the one-dimensional projections, supplied by the measurement 
system used, of the two-dimensional distribution of m(x, y) from numerous 
different angles of view. 

[00041 CT images can be generated both by means -wav of a CT unit with a 
scanning system that can circulate endlessly, and by meana -wav of a C-arc unit 
with a scanning system that can rotate only by less than 360°. The abbreviation 
"CT", for example in "CT raw data", is used below with reference to both types of 



unit 
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L00051 T he projection data are detemnned from the intensity I of an X-ray after its 
path through the slice to be imaged, and by its original intensity I 0 at the X-ray 
source in accordance with the absorption law 

In— ■ Jp(x,yJdI (1) 

1*0006] The integration path L represents the track of the X-ray considered through 
the two-dimensional attenuation distribution fx(x, y). An image projection is then 
compiled from the measured values, obtained with the aid of the X-rays from one 
direction of view, of the line integrals though the object slice. 

f00071 The projections originating from the most varied directions - characterized 
by the projection angle a - are obtained by a combined X-ray tube/detector system 
that rotates in the slice plane about the object. The units currently most employed 
are so-called "fan beam units" in the case of which tubes and an array of detectors 
(a linear an-angement of detectors) rotate together in the slice plane about a center 
of rotation, which is also the center of the circular measuring field. The "parallel 
beam units", which are affected by very long measuring times, will not be 
explained here, However, it may be pointed out that it is possible to convert from 
fan to parallel projections and vice versar^d^^jQiu^ m embodiment of the 
present invention, which is to be explained with the aid of a fan beam unit, can also 
be applied without restriction to parallel beam units. 

[00081 In fan beam geometry, a CT photograph c om prises includes line integral 
measured values -ln(I/I 0 ) of incident beams thai are characterized by a two- 
dimensional coupling of the projection angle as [0, In) and the fan angles 06 t^o, 
Pol defining the detector positions (p 0 being half the fan opening angle). Since rhe 
measuring system has only a finite number k of detector elements, and a 
measurement consists of a finite number y of projections, this coupling is discrete 
and can be represented by a matrix; 
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p(<V Pj : [o.2n] x [- p 0/ p 0 ) (2) 



or 



P(y,k):(l l 2...N P )x(l,2.„N s ) (3) 

The matrix p (y, k) is called a sinogram for fan beam geometry. The number of 
projections y and the number of channels k are of the order of magnitude of 1000. 

[00091 Taking the logarithms in accordance with the equation (1) thus yields the 
line integrals of all the projections 

p(a; 0) = In — = J/**, y) dl (3a), 

the totality of which is also denoted as radon transform of the distribution p(x, y). 
Such a radon transformation is reversible^^t is_thereby being possible to 
calculate m(x. y) from p(oc, 0) by back transformation (inverse radon 
transformation). 

lOOlOLIt is customary in this back transformation to apply a convolution algorithm 
in which the line integrals per projection are firstly convolved with the aid of a 
special function and then backprojected onto the image plane along the original 
beam directions. This special function, by m*mr-o&~ which the convolution 
algorithm is essentially characterized, is denoted as "convolution core". Owing to 
the mathematical configuration of the convolution core, the possibility exists of 
specifically influencing the image quality in the reconstruction of a CT image from 
the CT raw data. 

fOOll] For example, an appropriate convolution core can be used to emphasize 
high frequencies in order to increase the spatial resolution in the image. 
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Alternatively , or-a convolution core of an appropriately different type can be used 
to damp high frequencies in order to reduce the image noise. Thus, it may be stared 
in summary that when reconstructing images in eemptttefcomouted tomography 
the image characteristic, which is characterized by image sharpness/image noise 
and image contrast (two behaving in a complementary fashion relative to one 
another), can be influenced by selecting a suitable convolution core. There is a 
direct proportionality between image sharpness and image noise in this case, that is 
co say when the image sharpness is increased the noise is increased to the same 
extent. 

[0012] N o further investigation is now to be made into the principle of image 
reconstruction in CT by calculating the /i value distribution. A detailed description 
of CT image reconstruction is to be found, for example, in "Bildgebende Systems 
fur die medizinische Diagnostik" ["Imaging systems for medical diagnostics"], 3rd 
Edition, Munich: Publicis MCD Verlag, 1995, Publisher: Momeburg Heinz, ISBN 
3-89578-002-2 , the entire contents of which are hereby incorporated herein by 
reference . 

f0013T H owever, the task of image reconstruction is not yet concluded with the 
calculation of the // value distribution of the transiiradiated slice. The distribution 
of the attenuation coefficient fx represents in the field of medical application only 
an anatomical structure that still needs to be represented in the form of an X-ray 
image, 

f00141 In line with a proposal from G.N. Hounsfield, it has become generally 
customary to convert the values of the linear attenuation coefficient /i (which has 
the unit of measurement cm" 1 ) to a dimensionless scale in which water takes the 
value 0 and air the value -1000. The conversion formula to this "CT number" is: 

CT-number= ^—^^1000 (4) 
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rOQ!5l The unit of the CT number is called the "Hounsfield unit" (HU). This scale, 
denoted as the "Houndsfield scale", is very well suited to representing anatomical 
tissue, since the unit HU expresses the deviation in parts per thousand from 
and the fi values of most bodily substunces differ only slightly from the /x value of 
water. Only whole numbers are used from the number range (from -1000 for air up 
to approximately 3000) as carriers of image information. 

100161 However, the representation of the entire scale range of approximately 
4000 values would far exceed the discrimination threshold of the human eye. In 
addition, the viewer is frequently interested only in a small section of the 
attenuation value range, for example the differentiation of gray and white brain 
Substance, which differ only by about 10 HU. 

f00l71 Use is made of so-called image windowing for this reason. Only part of the 
CT value scale is selected in this case, and is spread over all available gray levels, 
Even small attenuation differences within the chosen window thus become 
perceptible gray tone differences, and all the CT values below the window are 
represented black and all the CT values above the window white. The image 
window can be varied arbitrarily both in terms of its central level and in terms of 
its width. 



f00181 _It is of interest in o o mputo r computed tomography to undertake multiplanar 
reformattings (MPR or secondary sections), Multiplanar reformattings are 
arbitrarily inclined, flat CT images calculated from a volumetric data record (also 
denoted as primary data record and usually represented by thin axial layers). Since 
the pixels generally do not occupy the position defined in the volumetric data 
record, and the layer thickness of an MPR is intended to be able to be set 
arbitrarily, it is necessary in this case to interpolate suitably. Particularly in the case 
of more recent CT units, the resolution of a volumetric data record is virtually 
isotropic. For this reason, it is possible to calculate from such a volumetric data 
record high-quality MPRs whose quality does not differ from that of the primary 
images. 
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£00191_However, it is likewise of interest during the course of a good diagnostic 
image evaluation to manipulate the image characteristic of an MPR - essentially 
characterized by sharpness and noise - by greens -use o f suitable filtere. Sharpness 
and noise of the reformattings are substantially determined by sharpness and noise 
of the primary, axial images as well as by the layer thickness set when generating 
the MPR. 

£0020]_The representation of CT images, in particular those of MPRs, with a 
different image characteristic is therefore of interest because a different evaluation 
(that is to say a different clinical assessment) of the same photograph of the 
corresponding th&us requires a different sort of representation of the recorded 
tissue. 

fO02]J_In the prior art, a targeted manipulation of the image characteristic of 
secondary sections is achieved by determining from the raw data a new volumetric 
data record by means-way of a new image reconstruction with the aid of changed 
convolution core parameters, and using this new primary data record to 
subsequently regenerate the original secondary sections. This means -refers to a 
reconstruction with the aid of a different convolution core that has a different 
characteristic, for example precisely a different sharpness. 

SUMMARY OF THE INVENTION 

r0022J_The above-described goal of adapting the image characteristic of CT 
images, in particular of secondary sections, optimally in accordance with the 
diagnostic requirements, can be realized only partly by mews-way of this known 
mode of procedure. A renewed first reconstruction with the aid of the parameters 
adapted to the later use changes the sharpness only transverse to the patient's axis. 
The patient's axis is denoted below as the z-axis. In an MPR generated from this 
reconstruction, the image characterizing properties remain unchanged in the z- 
direction. 
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f0023T In order also to achieve the desired success in the z-direction, the slice 
sensitivity profiles (SSPs) of this first reconstruction would have to be individually 
adapted. To date, this adaptation has been possible only in discrete seeps, that is to 
say a function for arbitrarily accurate tuning of the SSPs - which in principle can 
also be denoted as axial sharpness - is not provided in current CT units. Even if an 
arbitrary SSP could be used in the first reconstruction, the generation of suitably 
changed secondary sections would require a substantial outlay on operational and 
computer time. 

f00241 I t is therefore the^-an object of an embodiment of the present invention to 
propose* with regard to diagnostically optimized points of view A new techniques 
for improving and/or simplifying the filtering of CT images, 

f00251 T hus, an embodiment of the invention proposes a first inventive method for 
filtering CT images, having the following steps; 

a) acquiring a CT raw data record with the aid of a CT unit or of a C-arm unit, 

b) reconstructing a primary data record from the CT raw data record by means of, 
for example, a sharp convolution core and, for example, a narrow slice 
sensitivity profile, 

c) providing a transfer function as functional relationship between window width 
and image sharpness, and 

d) automatically calculating the image sharpness of the CT image of a selected 
slice, situated in the primary data record, as a function of a selected window 
width for the selected slice by means of an image processing procedure on 
the basis of the transfer function - 

fQ0261In this case, in accordance with a first embodiment of this first inventive 
method, in the image processing procedure the window width is coupled to at least 
one parameter of the three-dimensional convolution core with the aid of which the 
primary data record is reconvolved, and the at least one slice is modified with 
regard to the image sharpness. 
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f00271 I n accordance with a second exemplary embodiment of the first inventive 
method, in the image processing procedure the window width is coupled to at least 
one parameter of a two-dimensional convolution core with the aid of which the at 
least one slice is convolved and modified with regard to the image sharpness. 

f00281 In both exemplary embodiments, a renewed reconstruction with changed 
core parameters, and the renewed calculation of secondary sections are avoided^ 
Aand-4s such A this constitutes a substantially quicker method 

fOQ291 It is advantageous initially that the primary data record is reconstructed 
with the aid of a convolution core of maximum sharpness and a slice sensitivity 
profile of maximum narrowness. 

r00301 According to an embodiment of the invention, the slice can be an axial slice 
or else a secondary slice. 

f0_031] T he selection of the window width is performed in this case by the user by 
means of a mouse or keyboard. 

[0032LFurthermore, it can be advantageous when the transfer function additionally 
constitutes a functional relationship between window width and slice sensitivity 
profile. 

[003JL According to an embodiment of the invention, a second inventive method is 
further proposed for filtering CT images, having the following steps: 

(1) acquiring a CT raw record with the aid of a CT unit or of a C-arm unit, 

(2) reconstructing a primary data record with, for example, a sharp convolution 
core and, for example, a narrow slice sensitivity profile, 

(3) reconstructing an image stack with the aid of a corresponding image 
characteristic on the basis of the primary data record, 

8 
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(4) calculating a changed image characteristic of ihe image stack by means of 
an image processing procedure running in the background on the image 
computer, and 

(5) visualizing the image stack in the form of CT images with the aid of the 
changed image characteristic, 

[0034] In a, once agam, first exemplary embodiment of this second method 
according to the invention, via an input interface^ the user can vary at least one 
parameter of the three-dimensional primary data convolution core a which is then 
used to re^convolve the primary data record in the image processing procedurer-at*d 
■ F from this* a new image stack is determined that i$ modified with regard to 
image sharpness. 

f0_0351 In a second exemplary embodiment of the second method according to the 
invention, via an input interface^ the user can vary at least one parameter of a two- 
dimensional convolution cor e,_Thi s-wbisk is then used in the image processing 
procedure for the individual convolution of individual slices of the image stack* 
and to modify them with regard to their image sharpness. 

f00361 Steps (4) and (5) can be repeated by the user until a satisfactory image 
characteristic is achieved. 

I2037L According to an embodiment of the invention, the image stack can 
constitute a stack of axial images or else a stack of arbitrary secondary sections. 

f0038] A further aspect of an embodiment of the present invention is that the slice 
thickness of a slice contained in the image stuck can be varied in the event of an 
unsatisfactory image characteristic. 

JjM)39]JFurthermore, there is proposed^ according to an embodiment of the 
invention, a oomputor epmputed tomography unit for carrying out the method 
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which has a computer that carries out the individual steps of the CT image 
processing, 

fOQ401 L ikewise, there is proposed^ according to an embodiment of the invention 
a computer software product that implements a method when it runs on a 
computing device connected to a comput e r computed tomography unit, 

BRIEF DESCRIPTION OF THE DRAWINGS 

f00411 F urther features, properties and advantages of the present invention will 
now be explained with the aid of exemplary embodiments and with reference to the 
accompanying figures of rhe drawings, in which: 

figure 1 shows a schematic of a CT apparatus for a fan bea!m method in 
accordance with an embodiment of the present invention, 

figure 2 shows a Hounsfield scale in which the Hounsfield units (HE) of 
different organs of the human body are specified, 

figure 3 shows windowing in the case of the representation of CT images, 

figure 4 shows an inventive transfer function, 

figure 5 shows the equipotential surface of a 3D convolution core in the form of 
an ellipsoid that is cut by a characteristic plane, and 

figure 6 show 5 , in the form of an ellipse, the equipotential line given by the 
characteristic plane, 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
I0042LUlustrated in figure 1 4 schematically, is a comtmto r commited tomography 
unit for a fan beam method that operates in accordance with an embodiment of the 
present invention. In this device, X-ray tubes 1 and radiation detectors 2 rotate 
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jointly about a center of rotation that is also the center of the circular measuring 
field 5 and at which the patient 3 to be examined is located on a patient couch 4. 
Instead of the comp t rte r compmed tomography unit, it will aUo be possible to use a 
C-arm unit.-_The patient couch can be displaced along the longitudinal body axis 
in order to be able to examine different parallel planes of the patient 3. 

f0p_431 As may be seen from the drawing, CT photographs yield transverse 
sectional images, that is to say images of body slices that are oriented substantially 
perpendicular to the body axis. This method of representing slices constitutes the 
distribution of the attenuation value // s (x, y) itself (z is the position on the 
longitudinal body axis). Computer Computed tomography (termed CT below) 
requires projections at very many angles a In order to produce a tomogram, the 
radiation cone emitted by the X-ray tube 1 is masked out such that a plane beam 
fan results that cast one-dimensional central projections of the transirradiated slice, 

f0044T In order to reconstruct exactly the distribution of the attenuation values 
y)< this beam fan fR^e^ -should be perpendicular to the rotation axis^, and ? 
{Furthermore, it should be so widely spread that it completely covers the targeted 
slice of the measurement object from any projection direction a This beam fan 
penetrating the object is picked up by detectors that are arranged Imearly on a 
circular segment. There are up to 1000 detectors in commercially available units. 
The individual detector reacts to the incident beams with electric signals whose 
amplitude is proportional to the intensity of these beams, 

[0045] Each individual detector signal belonging to a projection a is picked up 
respectively by an electronic measuring system 7 and passed on to a computer 8. 
The computer 8 is used to process the measured data in a suitable way and firstly 
to visualize them in the form of a sinogram (in which the projection cx is plotted as 
a function of the measured values of the corresponding channel (3) in so-called 
Gordon units, but finally in the form of a natural X-ray image in Hounsfield units 
on a monitor 6. 
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100461 The-An aim of an embodiment of ihe present invention is to provide simple 
and quick methods in order to be able to adapt the image characteristic of CT 
images - in particular of secondary sections - from diagnostic points of view. 



[0047]_A first exemplary embodiment of a first inventive method con ot Ota r esides - 
by contrast with the above-described prior art - in undertaking direct filtering of the 
volumetric data record obtained from the raw data. The raw data are measured in a 
first step by mews-way of detector elements 2 in the form of attenuation profiles of 
different projectionSr-a«d-t^_These in turn are measured in different neighboring 
slices along the body axis. The raw data are picked up by the electronic measuring 
system 7 and passed onto the computer 8, which calculates from this CT raw data 
record in a second step 4 a primary data record (volumetric data record) that 
constitutes^ as such in the form of HU units* the basis for the calculation of axial or 
secondary sectional images. All further steps of the method according to an 
embodiment of.the invention are carried out in a purely computational fashion in 
the computer 8 of the CT unit anchor visualized on the monitor 6. 

£004|]_Thus, the inventive filtering of the primary data record is performed 
retrospectively in the first exemplary embodiment of the first inventive method by 
*»0a»s-way_of a three-dimensional convolution by applying a three-dimensional 
convolution core in the form of a three-dimensional matrix directly to the 
volumetric data record. The volumetric data record is thus changed in such a way 
that the secondary sections produced therefrom have the desired image 
characteristic. 



r<H>491 T he frequency response of the axial primary images in the transverse 
direction (that is to say in the x/y direction) is determined by the modulation 
transfer function of the reconstruction algorithm used and is generally isotropic. 
The modulation transfer function indicates which (spatial) frequency is visible in 
the image with which amplitude. The modulation transfer function is dependent 
essentially on the geometry of the X-ray system (focal length, number of the 
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detector channels in the transverse direction, number of the projections per cycle 
etc.), and of the convolution core that is used in the reconstruction. The SSP, 



detector elements etc.) and algorithm components (axial weighting function in the 
spiral algorithm), defines in the z-direction the behavior in the frequency domain. 
As a rule, this direction is distinguished from the axial slice planes - even in the 
case of the 3D isotropic resolution. 

r00501 A three-dimensional convolution core with the desired effect on the image 
characteristic can be radially symmetric with reference to the z-axis, while the z- 
axis assumes a special role, In general, the filter characteristic differs, however, 
along aJl three space axes. 

[00511 It may be assumed for the purpose of illustration that the 3D convolution 
core corresponding to the filtering is a Gaussian function, the width thereof along 
all the axes generally being different. The equipotential surface of such a 
convolution core is ihen the surface of an ellipsoid with different semi-axes, as is 
illustrated in figure 5, by way of example. Such an ellipsoid is used to determine 
the volumetric data record on the basis of which the secondary sections (MPRs) are 
then interpolated. 

[0052] A second exemplary embodiment of the first inventive method conoioto 
resides in directly filtering the secondary sections (MPRs). The usual procedure in 
this case is to undertake a two-dimensional convolution on the basis of the 
previously set slice thickness, as well as the spatial position of a secondary section. 
The two-dimensional convolution core required for this purpose is a two- 
dimensional matrix that is provided by the intersection of the 3D convolution core 
with the characteristic plane of the corresponding secondary section, The 
characteristic plane is uniquely fixed by a plane, parallel to the secondary sections, 
through the origin, and defines the position Of the secondary sections. In general, 
the result for the equipotential lines of the two-dimensional convolution core is 




including measuring system components (focal length, width of the 
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ellipses with axes that are rotated in general with reference to the Cartesian 
geometry. Such an ellipse is illustrated in figure 6. 

r005_33 Basically, both exemplary embodiments - 3D convolution of the volumetric 
data record and subsequent calculation of the MPRs (3D filter + MPR) and 
calculation of the MPRs and their subsequent filtering with the aid of a 2D 
convolution core (MPR + 2D filter) - are equivalent, to the extent that the 
properties perpendicular to the MPR plane are determined by the algorithm for 
calculating the secondary sections. Thus, for example, the slice thickness is set by 
the user during the calculation of an MPR with the aid of an interpolation function. 
Particularly in the case of two-dimensional filtering, the aim is also to change only 
the transverse characteristics of the MPRs. 

r00541 T he second exemplary embodiment, in particular, of the first inventive 
method (two-dimensional convolution core applied to MPRs already present) 
permits the image characteristic (sharpness, noise) of already present MPRs to be 
changed in a single step. In the case of typical computing times of approximately 
100 msec per MPR, this method is substantially faster than the conventional 
method, in which a renewed reconstruction of the primary data record and a 
subsequent renewed calculation of the secondary sections are carried out from the 
raw data. On the basis of a stack of 100 primary images, from which 50 secondary 
sections have been calculated, the conventional method requires 50 seconds purely 
for the reconstruction, assuming two reconstructed primary images per second. In 
addition, there is also the production of the secondary sections. Direct filtering in 
accordance with the second exemplary embodiment of the first method of the 
present invention requires, by contrast, only about 5 seconds, 

0)055] The exemplary embodiment of the first inventive method in the form of the 
3D convolution is not quite so fast as the direct convolution of the already 
calculated secondary sections, since the volumetric data record is predetermined 
with the aid of a 3D convolution core before a renewed calculation of secondary 
sections, and this is more time consuming than a pure convolution of existing 
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MPRs. The advantage of such a mode of procedure consists, however, in that 
arbitrary secondary sections that, in the final analysis, have the desired image 
characteristic can be laid in retrospect through the newly calculated volume. 

[OfiSgLA further aspect of the present invention is to provide the user in a second 
inventive method with a possibility of being able to adapt the image characteristic 
of a CT image (in particular an MPR) to a diagnostic requirement in a simple and 
quick way by «*da«e-use_of a suitable visual interface which is ultimately to be 
implemented in the computer 8. Specifically, it is thus possible to represent 
different anatomical tissue types with a different image characteristic (image 
sharpness and image noise). 

fOQS71 The second inventive method of an embodiment p f the present invention 
capitalizes on the fact that the CT values (HU values) occupy different ranges in 
the Hounsfield scale depending on the organic structure. 

[0058] The Hounsfield scale is illustrated in figure 2. The CT values of the 
individual organs assume specific ranges largely independently of the X-ray 
spectrum used. Thus, because of their low density and the Jow attenuation caused 
thereby, lung tissue and fat exhibit negative CT values in the range from *950 to - 
550 and -100 to -80, respectively. Most other tissues are in the positive range 
(kidney: 20-40, heart: 40-100. blood; 50-60, liver: 50-70). Because of the high 
atomic number of calcium and thus the relatively high attenuation, bone tissue has 
CT values of up to 2000 HU. 

£00591,4096 (= 2 12 ) different gray scale values are available for conventional CT 
units for the purpose of representing the overall Hounsfield scale. However, the 
viewer can distinguish at most 60 to 80 gray levels. For this reason - as already 
mentioned - in CT imaging use is made of windowing in which an HU-interva! of 
interest is assigned the entire gray scale. This is illustrated by way of example in 
figure 3 for the range of compact bone substance. The user interactively fixes the 
window in terms of center (window position C) and width (window width W), for 
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example by means of mouse or rotary knob. In the example of figure 3, the center 
is at C = 2000 and the window width is at W = 400. The window is assigned 10 
gray levels between white and black. 

100601 In general, the windowing is performed soch that the organic structure of 
interest is optimally expressed. A classic application is a CT lung photograph 
(thorax), in the ca$e of which both the contrasting bronchial branches and the low- 
contrast soft tissue of the lung have diagnostic relevance. In order to illustrate the 
bronchial branches, it is expedient to select a large window, the advantage thereby 
being that the noise is optically reduced, since the noise exhibits a maximum 
amplitude that is increasingly suppressed the larger the windowing. In order to 
resolve the fine structure of the tissue of the soft part, it is suitable to select a small 
window in which although the image sharpness is high a high noise amplitude 
associated therewith is also accepted at the same time, however, 

f00611 It is sensible from a diagnostic point of view when the window width 
selected for representing the image fixes the sharpness parameters and thus the 
image characteristic. The inventive idea thus consists in coupling the window 
width selected for the image representation to image sharpness and noise. For this 
purpose, an embodiment of the invention defines a transfer function that uniquely 
assigns a convolution core to the selected window width. Such a transfer function 
is illustrated in figure 4. The abscissa represents the window width, and the 
ordinate the image sharpness. The transfer function is advantageously a 
monotonically rising curve, since with increasing window width for a given 
maximum noise amplitude the noise is suppressed more and more and the image 
contrast is increased. In general, however, it is possible to select between the 
window width and image sharpness an arbitrary relationship that expediently 
supports the diagnosis. 

f0062 XOwjng to this functional relationship between window width and image 
sharpness, which is implemented, according to an embodiment of th* invention, by 
software in the computer 8 of the CT unit, in accordance with a first exemplary 
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embodiment of the second inventive method the user is able to change the image 
sharpness of the CT image of a selected slice present in the volumetric data record 
by filtering this basic volumetric data record in an image processing procedure 
based on the transfer function, As already described in detail above, it is possible 
in the general case for the filter characteristic (the nature of the convolution core 
on which the filtering is based) to differ along all three space axes. Clearly, the 
cquipotential surfaces of the convolution core are then surfaces of ellipsoids with 
different semi-axes. Filtering the volumetric data record has the advantage that 
arbitrary slices can be obtained in retrospect in this volumetric data record with the 
correspondingly modified image characteristic. 



r0063] fo r reasons of (computational) speed, however, in a second exemplary 
embodiment of the second inventive method direct filtering of the secondary 
sections is undertaken by using a two-dimensional convolution core. The spatial 
position of the secondary sections (MPRs) is uniquely fixed by a plane, parallel to 
the secondary sections, through the origin of coordinates (denoted below as 
"characteristic plane"). As already explained and illustrated in figure 6, the 2D 
convolution core used is given as a section of the 3D convolution core (ellipsoid) 
with the characteristic plane in the forai of an ellipse. 

fOQ641 The filtering according to the second exemplary embodiment of the second 
inventive method is implemented by carrying out a two-dimensional convolution 
on the corresponding secondary plane (MPR) of a prescribed reference volume. 
The reference volume is the primary data record (original volumetric data record) 
that was calculated from the raw data record, for example, with the aid of a 
convolution of maximum sharpness as well as a slice sensitivity profile of 
maximum narrowness. 



[OggS] In general, filtering of the above type according to the first exemplary 
embodiment of the second inventive merhod is carried out such that after 
prescription of the transfer function the three-dimensional filter characteristic is 
adapted as a function of the window width by means of an image processing 
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procedure that is carried out on rhe computer 8. The reference volume is adapted 
with the aid of the corresponding three-dimensional convolution core to the desired 
image characteristic, and the corresponding secondary planes (MPRs) are 
subsequently calculated therefrom, Finally, the respective secondary plane is 
displayed with the aid of the selected window width. 

[00661 The practical implementation of such window-controlled filtering of CT 
images is performed, for example, by virtue of the fact that the user changes the 
window width (window setting) in accordance with the diagnostic aims on the 
monitor 6 of the computer 8 at a virtual interface with the aid of the mouse, and 
correspondingly simultaneously varies the filter characteristic of the basic 
convolution core in accordance with the transfer function. Thus, it is possible 
solely by varying the window width to adapt the image characteristic (specifically, 
sharpness and noise) to the diagnostic requirements of the user. 

TO0671 Both inventive methods described are important not only for the adapted 
representation of secondary sectional images, but also in clinical practice, 
particularly for the diagnostic evaluation of axial slice images. 

100681 The invention being thus described, it will be obvious that the same mav be 
varied in many wavs. Such variations are not to be regarded as a departure from 
the spirit and scope of the inv ention, and all such modifications as would ft e 
obvious to one skilled in the art are intended to be included within the scope of the 
following claims. 
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Patent Cloimo WHAT IS CLAIMED IS: 

1. A method for window-controlled filtering of CT images, having - tho 
following stops comprising ; 

_a) enquiring a CT raw data record w&» the oid-of a CT unit or of a C orm unit, 

b)— _reconstructing a primary data record from *e-a_CT raw data record-by 
moans of for e xamples; q sharp convolution ooro and, for exampl e , a sorrow 
slice sensitivity profile^ 

, — providing a transfer function as functional relationship between 

window width and image sharpness-r^and 

&) — automatically calculating tbe-image sharpness of the CT image of at 

lease one selected image slice, situated in the primary data record, as a function 
of a selected window width for the selected at least one slice by moans o f usin^ 
an image processing procedure based on the baoiG of the transfer function. 

<# 

2. The method as claimed in claim 1, e horaotorfcod in that wherein. in the 
image processing procedore A the window width is coupled to at least one parameter 
of the-ajhree-dimensional convolution core* with the aid of which the primary data 
record is reconvolved, and wherein the at least one selected slice is modified with 
regard to the image sharpness. 

3. The method as claimed in claim 1, choniot o rigod in thci t vvheran. in the 
image processing procedure^ the window width is coupled to at least one parameter 
of a two-dimensional convolution core* with the aid of which the at least one slice 
is convolved and modified with regard to the image sharpness. 

4. The method as claimed in claim I, chDraoteFK e44fl^ha twherein the primary 
data record is reconstructed with the aid of a convolution core of maximum 
sharpness and a slice sensitivity profile of maximum narrowness. 
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5. The method as claimed in claims 1*^4, characterised in tha t wherein the art 
least one slice is an axial slice. 

6. The method as claimed in claims characterised in tha t wherejfl the at 
least one slice is a secondary slice. 

7. The method as claimed in claimo 1 to 6, ohoraotoriaod in that claim 1, 
wherein the selection of the window width is performed by ihe-^useri by moana of 
ft- mous e or keyboard. 

8. The method as claimed in claim L wherein claims 1 to 7. ch a racterized i n 
that the transfer function additionally constitutes a functional relationship between 
window width and slice sensitivity profile. 

9. A method for retrospective filtering of CT images, ha ving the following 
step scomprising : 

W acquiring a CT raw data rooord wi t h th e aid of a CT unit or of a C 

are unit, 

(3) reconstructing a primary data record from a from a CT raw data record: eft 

the baoiG of, for example, a sharp convolution coro and, for examplo, a 

narrow slice oonoitivity profil e , - 
t 3 ) reconstructing an image stack with the aid of a corresponding image 

characteristic on the basis of the primary data record^ 
calculating a changed image characteristic of the image stac k by moans 

of an image processing procedure running in the background on the imago 

oomputcft -Land 

visualizing the image stack in the form of CT images with the aid of 

the changed image characteristic. 

10. The method as claimed in claim 9, characterised in thatw herein via an input 
interf a ce tho user can vary at least one parameter of the three-dimensional primary 
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data convolution core is varied, w hteh-thcis then uood to rooonvolvo the- primary 
data record is then reconvolved in the image processing procedure using varied 
parameter , and fr o m thfc -a new image stack is determined from the reconvolved 
primary record thai is modified with regard to image sharpness. 

11. The method as claimed in claims 9 and 10 , charact e rized in that via an 
input interface the uoor con - vorv wherein at least one parameter of a two- 
dimensional convolution core is varied, the varied at least one parameter wfr»eh-is 
then used in the image processing procedure for the individual convolution of 
individual slices of the image stack and to modify them with regard to their image 
sharpness. 

12. The method as claimed in claims 9 to 11 , e h&raetorizod in tha t wherein the 
steps of calculating and visualizing ^) and (5) can bo ropoat e d by tho noar ^e 
repeatable until a satisfactory image characteristic is achieved, 

13. The method as claimed in olaimo 0 to 12, characteriz e d in thnt claim 9. 
wherein the image stack constitutes a stack of axial images, 

14. The method as claimed in claims 9 4o 13, characteris e d in tha t , wherein the 
image stack constitutes a stack of arbitrary secondary sections. 

15. The method as claimed in eloimo 9 to 11, characterized in tha t claim 9. 
therein *e-slicc thickness of a slice contained in the image stack is also varied in 
the event of an unsatisfactory image characteristic. 

16 - A o ompu t e r computed tomography unit for carrying out the method as 
claimed in one of claims 1 to IS cJaim 1 , t^vmg <omprising a compute finE device 
that oarrios out tho individual stops of th e CT imago processing . 



17. A computer software product e horooteriss e d - m that it implcmonts adapted to 
cause a computing device connected to a computed tomography unit to implement 
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the method as claimed in claim - n mothod a s claimed in ono of claims 1 to - 1 5 _ 1 
when it mno on a computing dovioo conn e cted to a oomputor tomography unit . 

18. The method as claimed in claim L farther comprising: 

acquiring the CT raw data record with the aid of at least one of a CT unit 

and a C-arm unit. 

19. The method as claimed in claim L wherein the primary data record is 
reconstructed from the CT raw data record using at least one of a sharp 
convolution core and a narrow slice sensitivity profile. 

20. The method as claimed in claim 18, wherein the primary data record is 
reconstructed from the CT raw data record using at least one of a sharp 
convolution core and a narrow slice sensitivity profile, 

21. The met hod as claimed in claim 2. wherein the at least one slice is an axial 

22. The meth od as claimed in claims 2. wherein the at least one slice is a 
secondary slice. 

23. The meth od as claimed in claim 3, wherein the at least one slice is an axial 

Slice. 

21 The method as claim ed in claims 3, wherein the at least one slice is a 

secondary slice. 

The method as claime d in claim 4, wherein the at least one slice is an axial 

slice. 

2& The method as claim ed in claims 4, wherein the at least one slice is a 

secondary slice. 
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27. The method as claimed in claim 1, wherein the selection of the window 
width is performed by a user using at least one of a mouse and keyboard. 

28. The method asclaimed in claim 9, further comprisin g: 

acquiring the CT raw data record with the aid of at least one of a CT unit 

and a C-arm unit. 

29. T he method as claimed_in claim 9, wherein the primary data, record is 
reconstructed from the CT raw data record using at least one of a sharp 
convolution core and a narrow slice sensitivity profile. 

30,, T he method as claimed in claim 28. wherein the Primary data record is 
reconstructed from t he CT raw data record using at least one of a sharp 
convolution core and a narrow slice sensitivity profile. 

1L The method as claimed in claim 9. wherein a changed image characteristic 

of the image stack is c alculated bv an image processing procedure running in the 
background on the i mage computer, 

32j The method as claim ed in claim 30, wherein a changed image characteristic 

of_the image stack is calculated bv an image processing procedure running in the 
background on an image computer. 

^ The method as claime d in claim 9. wherein, via an interface; at least one 

gara meter of the three-dimensional primary data convolution core is varied bv a 
usey, the primary data record is then rcconvolved in the image processin g 
procedure using varied parameter, and a new image stack is determined from the 
^convolved primary record that is modifi ed with regard to image sharpness. 

^ The rqet hod as claimed in claims 10 t wherein the steps of calculating and 

visualizing are repeatable until a sati sfactory image characteristic is achieved. 
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The method as claimed in claim 10. wherein the image stack constitutes a 

stack of axial images. 

36, Th e method as claimed in claims 10. wherein the image stack constitutes a 
stack of arbitrary secondary sections. 

37. T he method as claimed in claims U, wherein the steps of calculating and 
Visualizing are repeaiable until a satisfactory image characteristic is achieved. 

3S,_ The method as claimed in claim IL wherein the image stack constitutes a 
stack of axial images. 

39. The method a s claimed in claims 11. wherein the, image stack constitutes a 
stack of arbitrary secondary sections. 

40. A computed tomography unit for carrying out the method as claimed in 
claim 9. comprising a computing device. 

A compu ter software product adapted to cause a computing device 

connected to a comp uted tomography unit to_ implement the method as claimed in 
claim 9. 



4 2 . A computed tomography unit for carry ing out the method as claimed in 
claim 20, comprising a computing device. 

£L A computer software product adapted to cause a computing device 

connected to a computed tomogr aphy unit to implemencjhe method as claimed in 
claim 20. 



^ A computed tomography unit for car rying out the method as claimed in 

claim 32. comprising a computing device. 
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41 A computer software product a dapted to cause a computing dev ice 

connected to a computed tomography unit to impleme n t the method as claimed in 
claim 32. 
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ABSTRACT 

Re tronpootivo and/or window controlled filtering of imagoa for adapting the 
oharpnoofi - QTid noise in o o mputor tomography 

Th e present invention rotatoG to a 

A method is_fo r r fot oxamplo, window-controlled filtering of CT images , having 
t h eH ToHowinE otopfe g _ A 

dh- acquiring -^ CT raw data record is_acqui red, for example, with the aid of a 

CT unit or of a C-arc unit. T 

h) A _roconotni o ting a primary data record is then reconstructed f rom the CT 

raw data record by meansway of, for example, a sharp convolution core and, for 
example, a narrow slice sensitivity profi)e^A 7 

providing n transfer function is then provided as functional relationship 

between wmdow width and image sharpnes$7-em4, 

automatical l y calcul ati ng ■ tl he image sharpness of the CT image of a 

selected slice is then automatical ly calculated , situated in the primary data record, 
as a function of a selected window width for the selected slice by me*fts -wav of an 
image processing procedure on the basis of the transfer function. 

Figure 4 



